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Abstract 

An  engineering  model  of  a  steady  capillary  plasma  discharge  is  presented.  Two 
regimes  of  operation  are  studied  and  compared,  differing  by  the  assumptions  made  in 
the  treatment  of  the  radiative  fluxes.  In  the  first  case,  the  plasma  density  is  high 
enough  that  the  radiation  Rosseland  mean  free  path  is  smaller  than  the  capillary 
radius;  this  case  of  black-body  radiation  leads  to  super-high  pressures.  In  the  second 
case,  the  radiation  mean  free  path  is  larger  than  the  capillary  length,  and  the  limit  of 
volume  radiation  is  achieved,  leading  to  more  moderate  pressures.  Both  assumptions 
lead  to  different  scaling  laws,  with  consequences  for  the  capillary  design.  The  current 
model  is  based  on  the  conservation  laws  for  mass,  momentum  and  energy,  and 
assumes  local  thermodynamic  equilibrium  (LTE)  and  sonic  conditions  at  the  exit 
plane.  Various  assumptions  have  been  made  in  the  literature  about  the  amount  of  heat 
being  used  for  the  phase  change  and  the  ionization  cost  of  the  wall  material,  and  the 
effect  of  these  assumptions  is  also  examined.  Finally,  we  discuss  a  number  of 
potential  applications  of  the  discharge,  describe  planned  experiments,  and  discuss  on¬ 
going  work  in  extending  the  model  to  multi-dimensions  and  the  intermediate 
(between  black-body  and  volume  radiation)  regime. 

1  Introduction 

The  capillary  discharge  has  been  previously  studied  as  potential  electro-thermal  and 
electro-thermo-chemical  guns  device  by  various  investigators  since  the  mid  1980s, 
e.g.  Tidman  et  al.  (1986),  Loeb  and  Kaplan  (1989),  Powell  and  Zielinski  (1992),  Zoler 
et  al.  (1994,  1995),  Gilligan  and  Mohanti  (1990),  and  Shafir  et  al.  (2005).  These 
studies  investigated  the  dynamics  of  high-pressure  ablative  plasma  discharges 
achieved  through  the  dissipation  of  high  (several  kA)  current  pulses  in  a  narrow 
channel.  Other  types  of  capillary  discharges  operating  with  much  faster  pulses  and  at 
lower  pressure  have  more  recently  been  the  subject  of  increased  attention  as  potential 
devices  for  X-ray  lasers,  wakefield  acceleration  and  ultra-intense  laser  guiding  (e.g., 
Cros  et  al.,  2000).  Therefore,  the  capillary  discharge  can  operate  at  least  in  two 
regimes;  as  a  very  high-pressure,  ablative  plasma  discharge  and  a  low-pressure,  non¬ 
ablative  fast-pulse  discharge.  The  intermediate  regime  has  not  been  thoroughly 
explored,  but  is  also  of  potential  interest  for  propulsive  or  other  applications  with  a 
high  duty  cycle.  The  high-pressure,  low-temperature  regime  -  of  the  order  of  Giga- 
Pascal  (GPa)  and  1-3  eV  respectively  -  is  a  result  of  fast  ablation  of  the  capillary 
walls,  and  presents  formidable  challenges  for  re-usability,  due  to  the  very  high 
stresses.  The  low-pressure,  high  electron  temperature  -  of  the  order  of  Mega-Pascal 
(MPa)  and  100  eV  -  regime  is  more  amenable  to  repeated  use,  but  has  a  reduced 
thrust  or  acceleration  performance,  and  thennal  non-equilibrium  effects  can  be 
significant.  Ideally  the  intermediate  pressure  and  temperature  regime,  typically  of  the 
order  of  5-15  MPa  and  4-10  eV,  could  be  handled  with  current  material  strength,  and 
reduced  ablation  rates  could  allow  continuous  or  repeated  feedstock  injection 
(propellant).  This  approach  was  considered  by  Burton  et  al.  (1990,  1991)  for  the 
design  and  testing  of  a  prototype  Pulsed  Electrothennal  Thruster  (PET)  operating  at 
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30-50  MPa  and  2  eV,  yielding  a  thrust-to-power  ratio  77P  ~  lOmN/kW  with  water 
injection,  a  very  respectable  figure  of  merit  for  high-power  electric  propulsion  (EP) 
systems.  The  current  study  is  partially  motivated  by  the  need  to  further  explore  this 
intennediate  pressure  regime  and  the  need  to  better  evaluate  and  optimize  device 
performance  for  applications  of  interest.  For  example,  the  pressure  and  temperature 
(P,T)  conditions  mentioned  above  could  potentially  lead  to  a  thrust  density  of  at  least 
10  MN/nr  if  operated  continuously  and  a  specific  impulse  of  3,000  sec;  this 
corresponds  to  a  thrust  of  approximately  1  kN  per  capillary,  an  impressive  potential 
for  a  plasma  thruster  device.  Capillaries  can  be  combined  together,  improving  the 
thrust/mass  ratio,  but  the  actual  performance  would  probably  be  limited  by  the  duty 
cycle  than  can  be  practically  maintained.  Continuous  or  high-frequency  operation 
may  be  difficult,  due  to  the  high  rates  of  energy  deposition  and  heat  transfer  in  a  small 
volume;  therefore,  predictive  analysis  of  the  thennodynamics  of  the  discharge  is  an 
essential  first  step  in  the  development  of  this  device.  For  this  purpose,  we  have 
developed  a  zero-dimensional  (0D)  model  of  the  discharge  to  derive  scaling  laws  and 
global  properties  of  plasma  capillary  thrusters;  this  model  will  be  further  expanded  to 
one-  and  two-dimensional  models,  currently  under  development,  with  additional 
physics  and  unsteady  effects. 


2  Description  of  the  Model 

2,1  Conservation  Laws  and  Approximations 

The  ablative  capillary  discharge  can  be  approximately  described  by  a  three-layer 
configuration.  The  outer-most  layer  is  the  solid  wall,  usually  some  form  of 
polyethylene,  occasionally  Teflon  or  some  other  insulating  material,  which 
evaporates  and  enters  the  transition,  or  “boundary”  layer,  where  it  can  be  dissociated, 
ionized  and  heated.  The  innermost  layer  is  the  plasma  core,  usually  considered  fully 
ionized  and  uniform.  Electrodes  are  located  at  both  end,  with  one  being  open-ended, 
through  which  the  plasma  can  expand.  A  current  pulse  of  several  kA  -  from  2  kA  to 
140  kA  in  some  cases  -  and  of  milli-second  duration  -  from  0.1  to  5  ms  in  the 
ablative  case  -  flows  through  the  plasma  along  the  capillary.  The  large  Ohmic  heating 
is  responsible  for  rapid  temperature  rise,  ionization,  ablation,  and  radiation. 


The  first  zero-dimensional  analytical  model  of  ablative  capillary  discharges  was 
developed  by  Loeb  and  Kaplan  (1989);  using  a  number  of  simplifying  assumptions, 
these  authors  obtained  scaling  laws  for  the  plasma  temperature,  pressure,  wall 
ablation  rate,  and  other  parameters  of  the  capillary  discharge  as  function  of  the 
current,  radius  ( Rc )  and  length  (  Lc )  of  the  capillary,  a  significant  progress  the  time.  In 
later  works  -  see  for  example  Powell  and  Zielinski  (1992)  -  the  Saha  equation  was 
used  to  more  accurately  calculate  the  plasma  composition  and  enthalpy  of  the 
capillary  discharge.  The  most  important  assumptions  made  relate  to  the  amount  of 
radiation  emitted,  as  well  as  the  amount  used  for  ablation  and  phase  change.  In  a  0D 
model,  the  conservation  laws  for  mass  and  energy  can  be  expressed  as: 

dp 
dt 


■Pa  +  Pe 


d 

—  (/?£•)  =  77-7 
dt 


•  2 


+Pehe~  —  <lr 


Rc 


.  +  Pa 


(la) 

(lb) 


where  pa ,  pe  are  respectively  the  rates  of  mass  addition  by  ablation  and  removal  by 
flow  at  the  exit,  per  unit  volume.  Correspondingly,  the  mass  injection  and  removal 
( p a  >  0, p e  <  0 )  introduce  source  tenns  for  the  total  energy  in  equation  (lb).  Here,  we 
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use  the  following  notation:  the  specific  (per  unit  mass)  thermal  energy  and  enthalpy 
are  s  =  CvT/m,  h  =  C  T/m,  while  total  energy  and  enthalpy  include  kinetic  and 

internal  energies: 

s  =(/-l )_1  kT  tm  +ju:  +  A  sv  +AsD  +Asx  +A  2)  (2) 

and  similarly  for  the  total  enthalpy  (h  -s+kT/rn).  In  (2)  we  have  expressed  the 
specific  heat  at  constant  volume,  and  described  the  total  internal  energy  as  changes 
due  to  vaporization  (V),  dissociation  (D),  electronic  excitation  (X)  and  ionization  (I). 
If  the  plasma  was  comprised  on  molecular  species,  vibrational  and  rotational 
excitation  should  also  be  included;  however,  the  conditions  are  such  that  complete 
dissociation  can  be  assumed,  and  multi-stage  ionization  and  electronic  excitation 
make  the  bulk  of  the  internal,  or  phase-change  energy  A .  The  first  term  on  the 
RHS  of  equation  (lb)  is  the  Ohmic  heating,  and  is  determined  by  the  plasma 
conditions,  i.e.  the  resistivity  >/ .  The  third  term  is  the  radiation  flux  at  the  edge 
( r  =  Rc )  of  the  plasma  column.  Prior  capillary  models  have  used  various  assumptions 
to  relate  the  mass  and  enthalpy  flows  from  the  ablated  material  with  the  radiation  flux. 
For  example  in  the  time-dependent  OD-model  of  Gilligan  and  Mohanti  (1990),  the 
ablated  material  from  the  transition  layer  is  not  dissociated  or  ionized  before  entering 
into  the  plasma,  and  therefore  the  rate  of  ablation  is  detennined  by: 
paAsv  =(2!Rc)-qr\  ,  since  the  radiation  only  performs  the  work  of  vaporizing  the 

material.  Gilligan  and  Mohanti  also  assumed  the  thermal  radiation  qr  incident  on  the 
insulator  surface  to  be  a  fraction/ (the  grey  factor)  of  the  blackbody  radiation  leaving 
the  plasma,  qr=  f  ■  crSBT4  .  Assuming  local  thermodynamic  equilibrium  (LTE)  and 
using  the  Saha-Boltzmann  equation,  Gilligan  and  Mohanti  calculate  the  internal 
energy  of  the  plasma,  the  resistivity  and  other  thermodynamic  parameters.  To  get  an 
agreement  with  experiments  showing  very  small  smaller  ablation  rate,  they  must 
assume  that  only  a  small  fraction  ( /  « 1 )  of  the  radiation  is  absorbed  in  the  wall. 
This  model  differs  notably  from  that  of  Loeb  and  Kaplan  by  the  assumptions  made  in 
the  amount  of  radiant  energy  absorbed,  and  the  “work”  done  by  this  radiation.  This 
will  be  further  discussed  hereafter. 


The  0D  model  that  is  considered  here  is  similar  to  prior  work,  can  be  also  described 
by  equations  (la-b),  and  also  uses  some  key  approximations.  In  all  cases  considered 
here,  the  ratio  of  plasma  pressure  to  magnetic  pressure  p  is  always  large  and  MHD 
effects  can  be  neglected.  Viscous  effects  and  species  diffusion  can  also  be  neglected, 
due  to  the  high  flow  Reynolds  number,  typically  of  the  order  of  106.  Thennal 
equilibrium  is  obtained  between  electrons  and  heavy  particles,  i.e.  Te  «  Th ;  this  is 
characterized  by  the  LTE  factor,  the  ratio  of  the  thermal  energy  exchange  rate 
between  electrons  and  heavy  particles  to  the  electron  Ohmic  heating  rate,  which  for 
the  conditions  of  interest  is  typically  10  - 10 3 . 


2.2  Physical  Model 

The  electrical  conductivity  plays  a  key  role  in  the  model,  by  determining  the  rate  of 

Joule  heating.  Neglecting  the  induced  magnetic  field,  the  conductivity  is  given  by: 

2 


a  = 


n„e 


m  \v  +v  ) 

e  V ei  en  / 


(3) 


The  electron-neutral  collision  frequency  is: 
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(4) 


‘".I 


«TeJv. 
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with  ve  the  mean  electron  thermal  speed  and  <j^f  the  collision  cross-sections  for 
momentum  (m)  exchange  with  species  s;  these  are  taken  from  Powell  and  Zielinski: 


'  eC 


'  eH 


30  na 
11  n  a 


where  an-A7is„ti2lmae2^  5.29x10" 


2.64x10 


-19 


[nr] 


(5  a) 

1.49  x  10~19  [m2]  (5b) 

[m]  is  the  classical  Bohr  radius.  The  electron-ion 


collision  cross-section  is  defined  by  the  Spitzer  equation  (Spitzer  and  Harm,  1953): 

38 Z  n„e2 


v„,  —  - 


re  ■  meT 


3/2 


In  A 


(6) 


The  factor ye(Z)  in  equation  (6)  is  a  weak  function  of  Z  and  can  be  approximated 
as:  ye»0.58+0.1  ■  (Z  — 1) .  The  thermally-averaged  Coulomb  logarithm  In  A  can  be 
approximated  as:  InA  «  24.38  + 1.51n7jeV] -0.51nn  3  For  a  fully-ionized  gas,  the 

(Spitzer)  conductivity  can  be  written  as: 

/?  T3/2 

cx  =  4^L  with  /?*  1.07x10  2  (7) 

Z  In  A 

While  Loeb  and  Kaplan  introduce  a  temperature  and  density-dependent  (3  in  order  to 
use  (7)  for  non  fully-ionized  plasmas  as  well,  we  keep  equation  (3)  for  the  general 
case,  and  use  expression  (7)  for  the  fully-ionized  case  only,  which  will  be  useful  when 
considering  the  scaling  laws  for  a  simplified  model. 


The  plasma  will  be  considered  either  fully-ionized  or  in  chemical  equilibrium.  In  the 
latter  case,  the  electronic  partition  of  bound  levels  is  also  included.  We  consider  13, 
12,  and  11  pseudo-levels  -  groups  of  elementary  levels  with  combined  statistical 
degeneracy  and  averaged  energy  -  for  Ci,  Cn  and  Cm,  i.e.  neutral,  lst-ionized  and 
second-ionized  Carbon  atoms;  Hydrogen  (HO  is  modeled  by  6  pseudo-levels.  The 
Saha  equilibrium  function  is  given  by: 


Sm+l(T)  =  - 


n+ 1 


■  =  A-T 


3/2 


1  m+ 1 


-•e 


I  Cm 

kT 


(8) 


with  A  =  6.04  x  ] (P  [ m  'eV  3/2]  and  um  the  electronic  partition  function: 


u  =  g°  +g1 

m  o  m  o  m 


A E1(m 
kT 


A E„ 


+  g;n-e  k''  + .  (9) 

One  can  then  express  the  ion  densities  as  function  of  total  species  density  and  the 
equilibrium  constants.  For  Carbon, 

n>^) 

y>l  (10) 


tlm(C)  ~  nC 


(C) 


) 


where  m  indicates  the  ionization  stage.  These  expressions  are  similar  to  equation  (17) 
of  Powell  and  Zielinsky.  Equations  for  Hydrogen  are  obtained  similarly,  with  the 
exception  of  having  only  one  ionized  state.  Using  charge  neutrality  the  electron 
density  can  be  written  as: 


”e 

s  m(s ) 


ngwy> 


(11) 
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Solving  this  non-linear  equation  for  the  electron  density  as  function  of  the  temperature 
will  allow  us  to  determine  all  atom  and  ion  concentrations. 


The  plasma  is  also  considered  weakly  non-ideal,  with  the  main  effect  being  a 
lowering  of  the  ionization  potential  and  pressure  from  the  ideal  plasma  conditions. 
Thus,  in  the  model  by  Griem  (1962),  I  m  =  F^eal  -AI(n,T)  and  P  =  nkT-AP(n,T ) ,  with: 

f  V/2 
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me 
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4  nEoA.De 
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kT 


24  nk 
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1/2 
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[eV] 
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v1020[nT3] 
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where  XDe  is  the  Debye  length.  We  expect  the  lowering  to  be  as  large  as  A IQ  ~ 
AP  ~  2  MPa  for  the  discharge  conditions  considered  here,  i.e.  up  to  10%  effect. 


(12a) 


(12b) 


1.15  eV  and 


2.3  Model  Formulation 

We  consider  two  classes  of  discharges:  1)  a  high  pressure/optically  thick  case,  where 
the  average  (Rosseland)  radiation  m.f.p.  is  much  smaller  than  the  capillary  radius  and 
the  plasma  radiation  can  be  modeled  by  a  blackbody;  2)  a  low-pressure/optically  thin 
case,  where  the  radiation  m.f.p.  is  larger  than  the  characteristic  dimensions  and  the 
plasma  radiation  is  given  by  volumetric  radiative  emission.  These  assumptions  are 
limiting  cases  of  a  more  complex  radiation  transfer  model  required  for  a  more  detailed 
evaluation  of  the  radiative  flux  absorbed  by  the  transition  layer  and  the  capillary  wall; 
this  detailed  model  is  currently  under  development  and  will  be  described  in  the  future, 
while  the  current  approximate  model  provides  preliminary  estimates  of  the  heat 
transfer  regime  and  bounding  limits. 

The  first  version  of  our  0D  model  is  similar  to  Loeb  &  Kaplan.  In  this  model,  the 
plasma  radiation  ablates,  vaporizes,  ionizes  and  heats  the  ablative  gas  in  the  transition 
layer  up  to  plasma  conditions.  Also,  it  is  assumed  that  the  average  energy  cost  per  ion, 
that  includes  the  vaporization,  dissociation  and  ionization  energies,  is  a  fixed  constant, 
and  all  carbon  ions  are  twice  ionized.  According  to  this  model, paha  =  2qrad IRC,  and 
the  enthalpy  injection  and  radiative  losses  cancel  each  other  out  in  the  energy 
conservation  equation  (lb).  The  surface-integrated  fluxes  are  then  defined  for  the 
black-body  (BB)  and  volumetric  radiation  (VR)  cases: 

-  Optic  ally-thick  (BB)  case:  Fbb  =  (2ttRcL)  ■  crSBT4  (13a) 

-  Optic  ally-thin  (VR)  case:  FVR  =  (kR2cL)  ■  aVR  (Z  Z2 )  •  nfTul  (13b) 

where  in  SI  units,  <jvr  ~  1 .42  x  1 0  411 .  The  volumetric  radiation  in  (13b)  is  assumed  to  be 
due  to  continuous  emission  (bremsstrahlung),  assuming  that  all  line  radiations 
(including  bound-free  transitions)  are  self-trapped.  This  model  simplifies  the  scaling 
analysis,  while  the  free-bound  emission  will  be  added  later.  The  Z-factors  in  (13b) 
originate  from  a  factor  ^  nertkZ%  in  the  expression  of  the  bremsstrahlung  radiation, 

where  the  summation  runs  over  all  ion  species. 

Neglecting  the  kinetic  energy,  the  total  enthalpy  of  the  ablated  flow  is: 

K  ~  M~l  [a^  +Cp  (1  +Z){ac+aH  )  •  Vf]  (14) 
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Here,  ac  and  a  H  are  the  numbers  of  carbon  and  hydrogen  atoms  in  the  “polyethylene 
molecule”  of  mass  Ma  =  (acAc+aHAH)  -Mamu  (in  SI  units),  M amu  is  the  atomic  mass 
unit  and  C  =y l(y- 1)  is  the  specific  heat  at  constant-pressure,  withy  =  5/3.  As  in 


equation  (2),  A includes  the  energies  of  dissociation,  electronic  excitation  and 
ionization.  Using  the  fonnation  energies  of  the  atoms  and  molecules  in  various 
phases,  one  can  construct  the  chain  of  internal  energy  changes  for  each  ablating 
molecule  as  (for  simplicity  we  include  the  electronic  excitation  into  ionization): 


A  £0  = 


AfK]+  [A*, 


?CH(g)  SCH(s) 


b 

V 


A^x-/ . 


acsc  +  ocHoH 


CH(g) 


Jf  [ac (£°c++  - £°c)+aH (s°H+  -£°h)\  (15) 


—  acs c++  +  aH£R+  £ch(s)  ~  ac£c++  +  aHs H* 

Defining  an  average  formation  energy  of  ions  as: 

CLq£  r++  +  CC  [_f£  + 

« *  =  - ^LjL~ 

CLq  +  cc  h 


we  can  now  write  (14)  as: 

ha  «(«C+^)[M<I)+(1+Z)c^r] 

The  inflow  of  material  can  be  written  as: 


(16) 

(17) 


t-<  rad  r-r  rad 

V-pa=^-«Ma—  (18) 

ha  K 

with  Fmd  taking  one  or  the  other  form  described  in  (13a-b).  Defining  the  ion-acoustic 
speed  in  the  plasma  as: 

c(T)  =  ((1+Z)  •  ykTIM^f12  (19) 

where  Mt  =  AM amu  is  the  average  ion  mass,  the  outflow  through  the  open  end  of  the 
capillary  can  be  described  by: 

V-pe=-{nR1c)-clnl-Mi  (20) 

where  c,  =  c(T] ) ,  and  the  suffix  “1”  denotes  conditions  at  the  exit  plane,  where  sonic 
conditions  are  assumed.  In  the  quasi  steady-state,  the  rate  of  mass  injection  by 
ablation  is  approximately  equal  to  the  outflow  rate  (pa+pe~ 0).  Substituting  the 
plasma  radiation  fluxes  from  (13),  and  using  the  isentropic  flow  relations  we  obtain: 

(2.tzRcL)  -  {gsbt  )  =  -  {7iR2 ) ,  ^  .  r  ■  c(T0 )  (21a) 


-  (kR;L)-(cjvr(Z 


±  =  Mi(nRt)n0-r-c(T0)  (21b) 


where: 


f 


O+i) 


r  = 


2(y-l) 


\T+X, 

stands  for  stagnation  conditions, 


(22) 


a,  =njn,,  is  the  ion  fraction.  The 


the  suffix  “O’ 

volume-averaged  conditions  can  be  approximated  by  the  stagnation  conditions  in  the 
limit  of  a  long  capillary;  in  that  case,  the  suffix  “0”  is  dropped  and  replaced  by  “BB” 
and  “VR”  for  the  two  limiting  radiation  models  being  considered.  Assuming  a  fully 
ionized  plasma  (a,  =1),  equations  (21a-b)  can  then  be  inverted  to  yield  the  ion 
density  as  function  of  temperature: 
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i  2  L 

nBB=r - c 


f  T  V1/2 


nVR  _  r  '  ' 


_  c  (^) 
i?c  Cre/  t^  +  a+z)^^]  ^rre/; 

1  [Wo  +  (i+z)r  kT\  f  T  v/2 


(aVR(ZZ2)TV2) 


yT'ref  ) 


(23a) 

(23b) 


where  Tref  is  a  reference  temperature  (e.g.  1  eV)  and  cref  =c(Tref) .  In  the  energy 

equation  (2b),  and  since  it  was  assumed  that  paha  =2qrad  /Rc,  the  quasi  steady-state 
implies  that  the  Joule  heating  balances  the  enthalpy  outflow  at  the  open 
end ,(7rR2)clnlhi.  Since  the  flow  is  isentropic,  ht  is  also  the  stagnation  enthalpy  and 
therefore  the  total  heat  outflow  is: 


Fout  =K2)rc 


f  \1/2 

T 


ref 


\Jref  J 


■  11^  +  (\+z)c pkr] 


The  total  energy  deposition  by  Ohmic  heating  is: 


Q 


L  I2 


(24) 


(25) 


«  )  o’ 

with  /  the  current,  and  a  the  conductivity  (3).  The  energy  balance  (Q  =  FaEut  )  yields 
another  expression  for  the  average  ion  density  in  the  discharge: 

In  A 


n,  -  - 


T  T2  72 

L-l  ^  eff 


T-i  —1  —1 

t3/2_1  Cref 


f  T  ^ 


-1/2 


UQ+(\+Z)CpkTY 


(26) 

n1^  Z  pT'"  ™  1  L  ~  '  '  "  J 

This  expression  is  independent  of  the  assumption  made  about  the  Rosseland  radiation 
m.f.p.  and  can  be  applied  to  both  cases;  equality  of  (23a-b)  and  (26)  yields: 


I  bb 


3/2 


( ln2Z 


\  ^eff 


pa 


SB 


(  j,'y2  \ 

v!nAy 


(27a) 


/  Tc 

1  VR  ^  L ref 


f  j  ^ 


\Jref  ) 


Prt 1 


\  a VR  (Z  z  )  J 


f  T  ^ 


VlnAy 


u^+{\+Z)CpkT\  (27b) 


For  each  limiting  case  of  optical  thickness,  there  is  a  system  of  two  equations  (e.g., 
23a  and  27a)  with  two  unknowns  («,  and  T  ),  given  the  total  discharge  current,  which 
can  be  solved  iteratively.  Alternative  solutions  for  the  pressure  are  obtained  using  the 
assumed  equation  of  state  (EOS):  p  =  (1 +Z)njkT . 


3  Results 

3.1  Scaling  Laws 

Scaling  laws  can  be  extracted  from  the  formulas  (23)  and  (27);  by  also  ignoring  the 
u([>  term  in  these  equations,  it  is  possible  to  obtain  simple  relationships,  i.e.: 

-  BB  case:  T  oc  (Rc^,u) -/4/n  P  oc  (Z,A~32/11}7 14/11  n  oz  (L R26/n)  ■  Il0/n  (28a) 

-  VR  case:  T  oc  (L1,4Rlc/2)-  71/4  Poz  (L~V2R;1)  -IV2  n  oz  (L3/4)71/4  (28b) 

The  geometrical  scaling  laws  are  identified  by  the  terms  in  parentheses,  and  were 
easily  verified.  The  approximate  scaling  between  temperature,  current  and  pressure 
obtained  from  equations  (23)  and  (27)  could  also  be  verified:  in  the  BB  case,  and 
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assuming  a  constant  u<tl  ~  20  eV,  the  profiles  of  Figure  1  are  obtained.  At  high 
temperature,  «0  can  be  neglected  and  the  scaling  is  consistent  with  the  Loeb  & 
Kaplan  approximation  given  by  equations  (28a-b),  in  which  case  P  oc  r3-5.  At  low 
temperature  (T  «  2  cV),//,„  dominates  and  the  pressure  scaling  changes:  P  x  T4'5 .  For 
the  VR  case,  the  scaling  changes  from  P  cc  T,  I  <xT  at  low  temperature  to 
P ccT2 ,1  ecT2  at  high  temperature,  neglecting  the  Coulomb  logarithm  variation. 


Figure  1:  Scaling  of  pressure  and  current  (renormalized  by  the  Coulomb  logarithm) 
with  respect  to  temperature;  BB  and  VR  models,  fixed  u (], . 

3.2  Equilibrium  Model 

In  version  2  of  the  0D  model,  the  thermodynamic  equilibrium  formulation  is  used.  In 
that  case,  the  plasma  composition  and  the  phase  change  energy  are  computed  for 
each  value  of  the  temperature,  though  a  sequence  of  iterations.  The  mass  ablation  rate 
is  obtained  as  function  of  the  radiative  flux,  and  therefore  the  plasma  temperature; 
given  a  mass  ablation  rate  and  a  mass  outflow  rate,  the  plasma  density  itself  changes 
until  convergence  is  obtained.  Since  the  plasma  density  is  itself  a  function  of  the 
temperature,  only  the  plasma  temperature  is  required  to  obtain  the  final  plasma 
conditions,  for  either  model  of  radiative  flux.  The  phase  change  energy  can  be 
compared  to  the  constant  values  assumed  in  the  version  1  of  the  model;  the  impact  of 
the  exact  value  of  «0  is  noticeable  mostly  at  low  (<2  eV)  temperature.  The  pressure 
and  current  in  the  optically  thick  (BB)  and  thin  (VR)  regimes  could  then  be 
recomputed,  and  the  scaling  laws  (23)  and  (27)  could  be  recovered  in  the  limit  of 
high-temperature,  but  not  at  low  temperature.  This  is  expected,  since  the  simplified 
(version  1)  model  assumed  a  fully-ionized  plasma  with  twice-ionized  Carbon  atoms; 
calculations  indicate  that  below  4  eV,  this  assumption  starts  to  break  down. 

In  fact,  there  is  also  a  threshold  below  which  no  solution  can  be  found  in  the 
optic  ally-thin  regime.  This  is  best  observed  by  considering  equation  (21b):  at  a  given 
temperature,  and  for  the  thermodynamic  equilibrium  conditions,  one  cannot  generally 
assume  that  the  ionization  fraction  is  unity.  Therefore,  the  LHS  of  (21b)  -  which 
describes  the  incoming  energy  flux  -  varies  as /if,  while  the  RHS  -  which  describes 
the  outgoing  flow  -  varies  as  nh .  For  a  low-temperature  plasma,  the  ion  density  is 
determined  by  the  Saha  formula,  in  which  case  nf  oc  nh  and  both  energy  flows  are 
linear  functions  of  the  density.  For  a  fully-ionized  plasma,  the  Saha  equation  is 
saturated,  ni  ~  nh  and  the  point  of  equilibrium  would  be  obtained  at  the  point  of 
intersection  between  a  parabola  and  a  straight  line,  when  both  flows  are  plotted  as 
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function  ofnlt .  Depending  on  the  location  of  the  point  of  inflexion  for  the  incoming 
mass  flux,  there  may  not  be  an  intersection  between  the  two  curves.  This  can  be 
demonstrated  by  plotting  the  two  energy  flows  as  function  of  the  total  density,  for  two 
cases  of  temperature.  As  shown  in  Figure  2,  there  is  no  intersection  for  1  eV,  and 
therefore  no  possible  solution.  Figure  2  shows  the  switch  in  the  inflow  curves  between 
the  n2h  scaling  at  low  density  (which  favors  full  ionization)  and  the  nh  scaling  at  high 
density  (which  favors  partial  ionization).  There  is  a  shift  in  the  energy  outflow  profile, 
as  expected  (lines  with  no  symbols),  but  it  is  hardly  noticeable  on  this  scale. 


Figure  2:  Incoming  and  outgoing  total  enthalpy  flows  for  two  temperatures  and 
VR  model,  as  function  of  total  density. 


Figure  3:  Rosseland  mean-free-path  versus  (a)  electron  density  and  (b)  plasma  temperature. 

3.3  Results  and  Conclusions 

In  the  absence  of  a  full  spectral  radiation  transport  model,  the  validity  of  each 
radiation  model  considered  so  far  can  be  estimated  by  evaluating  the  Rosseland  mean- 
free-path  of  the  radiation,  XR .  As  shown  in  Figure  3  for  the  black-body  model,  the 
condition  of  an  optically-thick  plasma  ( AR  «  Rc ,  with  Rc  ~  1  mm  )  is  satisfied  only  at 
the  highest  values  of  density  and  temperature;  this  leads  to  extremely  high  pressures, 
of  the  order  of  10  GPa  (100,000  atm).  In  the  VR  case,  the  computed  plasma  density 
is  confined  within  a  much  narrower  range,  and  the  mean-free-path  behavior  is 
significantly  different.  In  that  case,  the  optically-thin  condition  (XR  >LC  »RC,  with 
Lc  ~  10  cm  )  is  satisfied  even  at  the  lowest  temperatures  of  interest.  The  true  physical 
condition  may  be  somewhere  in-between,  and  this  will  be  investigated  in  the  near¬ 
future  with  a  more  detailed  radiation  transport  model,  currently  under-development. 
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The  model  of  volumetric  radiation  considered  so  far  -  see  equation  (13b)  -  accounts 
only  for  the  free-free  emission.  Free-bound  emission  can  also  be  included  using  a 
well-known  approximation  (Zeldovich  and  Raizer,  2002).  The  modified  rate  of 
emission  now  reads: 


FVr  —  (x  RqL)  ■  <?VR  (Z  Z  2  )  •  nfT 


1/2 


1 


2  — 
T 


(29) 


As  shown  in  Figure  4,  for  the  same  temperature  range  the  enhanced  emission  given  by 
equation  (28)  leads  to  lower  plasma  densities,  as  well  as  lower  voltage  and  current.  It 
would  seem  a-priori  contradictory  to  obtain  lower  densities  when  the  radiative  flux  to 
the  capillary  walls  is  increased,  since  this  radiative  flux  must  be  balanced  by  the 
enthalpy  flux  of  the  ablated  products;  however,  one  must  take  into  account  that  as  the 
radiative  emission  is  enhanced,  the  steady-state  plasma  cannot  be  sustained  at  the 
given  temperature  with  the  similar  or  higher  density.  Thus,  the  equilibrium  is  reached 
at  lower  ion  densities,  leading  to  correspondingly  lower  values  of  pressure  and  Ohmic 
power.  The  impact  is  quite  significant;  according  to  Figure  4,  the  same  pressure  (TO7 
Pa,  or  100  atm)  can  be  obtained  for  the  same  power,  but  the  enhanced  radiation  case 
corresponds  to  a  temperature  of  8  eV,  compared  to  1.5  eV  for  the  free-free  emission 
only.  While  these  results  must  be  considered  with  extreme  caution  -  given  the  very 
approximate  model  used  so  far  -  they  indicate  that  an  optimal  regime  for  propulsive 
applications  would  be  found  in  the  (partially)  optically-thin  regime. 


Figure  4:  Impact  of  free-bound  emission  on  key  operating  conditions  for  same  range  of  temperature; 
free-free  model  is  described  by  equation  (13b),  while  equation  (29)  includes  also  free-bound  emission. 


The  thennodynamic  equilibrium  model  having  been  verified,  it  can  be  applied  to  the 
radiation  models  considered.  One  can  make  several  observations  from  the  results 
(including  free-bound  emission  for  the  VR  case),  shown  in  Figure  5: 

(1)  The  difference  in  pressure  (and  therefore  plasma  density)  between  the  BB  and  VR 
models  is  up  to  3  orders  of  magnitude;  clearly,  the  VR  case  leads  to  conditions 
much  more  amenable  to  longer  pulses,  repeated  use,  and  higher  duty  cycles. 

(2)  A  similar  observation  can  be  made  for  the  voltage  supply  (see  Figure  5-d)  and  the 
current,  both  decreasing  by  at  least  an  order  of  magnitude;  this  leads  to  reduced 
power  dissipation  and  more  compact  circuit  and  pulse  forming  network. 

(3)  The  deposited  (Ohmic)  power  is  also  considerably  reduced  in  the  VR  case, 
facilitating  thermal  management  for  high  duty  cycles. 

(4)  The  mass  flow  rate  is  typically  1  kg/s  f  the  BB  case,  leading  to  a  recession  rate  of 
approximately  1  mm/s  by  ablation;  in  the  VR  case,  the  required  mass  flow  is 
sufficiently  small  that  liquid  mass  injection  is  possible,  using  a  few  %  of  the 
capillary  wall  area  as  injector  area. 
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(5)  As  mentioned  in  the  previous  section,  the  BB  model  is  barely  valid  (according  to 
the  Rosseland  m.f.p.)  at  high  temperature  and  pressure,  while  the  VR  model 
appears  to  be  a  consistent  model  throughout  the  explored  conditions. 


Figure  5:  Plasma  conditions  for  BB  and  VR  (including  free-bound)  models,  Lc  =10cm  ,  Rc  =  2  mm , 
thermodynamic  equilibrium  computation. 
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Figure  6:  Left:  comparison  of  thrust  versus  speed  of  sound  at  sonic  point  for  the  three  radiation 
models,  for  single  capillary  of  radius  Rc  =  2  mm ,  assuming  steady-state;  Right:  propulsion  efficiency 
(left  axis)  and  thmst/power  ratio,  assuming  expansion  to  Mach  3. 

The  prospects  for  propulsive  applications  can  be  evaluated  from  Figure  6,  which  show 
the  potential  thrust  (assuming  steady-state)  that  can  be  achieved  by  a  single  capillary 
of  radius  2  mm.  Steady-state  operation,  or  at  least  high-duty  cycle,  being  more 
realistically  achieved  in  the  engine  volumetric  radiation  case,  one  can  presume  that  a 
thrust  of  approximately  100  N  per  capillary  would  be  obtained  with  a  specific  impulse 
(Isp)  in  excess  of  2000  s;  multiple  capillaries  could  be  bundled  onto  a  single  engine, 
and  lower  molecular  weight  propellants  can  be  used  to  increase  both  figures  of  merit. 
Thus,  the  capillary  discharge  has  the  potential  for  simultaneously  high-thrust  and 
high-Isp  propulsion,  provided  sufficient  power  is  available.  Also  shown  in  Figure  6, 
propulsive  efficiency  (kinetic  to  electric  power  ratio)  can  be  in  the  40-50%  range, 
comparable  with  a  state-of-the-art  EP  thrusters  (e.g.  Hall-effect  thruster  -  HET)  and 
thrust-to-power  ratio  ranges  between  30-100  mN/kW,  also  comparable  to  the  HET 
and  consistent  with  the  70  mN/kW  quoted  in  Section  1  for  the  PET.  Further  model 
development  and  analysis  is  required  before  more  definitive  conclusions  can  be 
drawn.  In  particular,  we  plan  to  accomplish  the  following  in  the  near-future; 

-  Implement  a  detailed  radiation  model  with  radiation  transport 

-  Examine  the  unsteady  equations  and  investigate  stability  conditions 

-  Include  heat  conduction  into  and  within  the  solid  walls 
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